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Abstract 

An asymptotic expansion for the grand partition function of an ideal Bose 
! gas is obtained for the canonical ensemble with arbitrary number of particles. 

It is shown that the expressions found are valid at all temperatures, including 
the critical region. A comparison of the asymptotic formulas for fluctuations 
of the Bose condensate with exact ones is carried out and their quantitative 
agreement is established. PACS: 05.30.Jp, 75.30.Ds, 75.70-i 
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S-J ■ Introduction 

on 

The statistical mechanics of an ideal Bose gas [1, 2] is one of the simplest while at 
the same time fundamental subject areas of theoretical physics. The most impressive 
^ ■ result of the theory is the remarkable phenomenon of Bose-Einstein condensation (BEC), 

which was in fact predicted by Einstein [3] the accumulation of an unlimited number of 
noninteracting (and this is also nontrivial) particles in their quantum ground state. We 
note that BEC turned out to be the first exactly solvable model of a phase transition 
and is a demonstration of the possibility of quantum behavior on macroscopic scales. 
Unfortunately, there are no substances in nature that satisfy the conditions of BEC in 
thermodynamic equilibrium. Nevertheless, the idea of the the Bose condensate has turned 
out to be very useful for understanding the essence of such remarkable physical phenomena 
as superfluidity and superconductivity and has been fruitful in various conceptual specu- 
lations in condensed matter physics and quantum field theory. 

By the end of the last century experimental technique had been perfected so much 
that it became possible to restrain a cloud of polarized atoms of hydrogen or alkali metals 
in a magnetic trap long enough to cool it to very low temperatures ( ~ 20 nK . . . ~ 1 /iK ) 
[4,6-8] as a result, a large number of atoms (~ 10 3 in the case of rubidium and ~ 10 5 for 
sodium) have been collected in the lowest-energy quantum state. Although it is understood 



that it is not in a state of thermodynamic equilibrium, this is commonly regarded as the 
experimental preparation of a Bose condensate. 

If one has in mind not a strict thermal equilibrium but a quasi-equilibrium, then it is 
entirely acceptable to consider a collective of some quasiparticles whose mutual interaction 
is also weak, as a rule, to be a suitable object for experimental investigation of BEC. As 
we know, the condition for thermal equilibrium here requires that the chemical potential 
vanish ( // = ), while the BEC regime corresponds to an increase of the chemical potential 
to the lowest energy level (// — > Sq ). In the case when Sq is small, an equilibrium gas of 
quasiparticles is found at the threshold of Bose condensation. However, in quasiequilibrium 
(e.g., under the influence of the external pump creating the quasiparticles) it is possible in 
principle to control the chemical potential, increasing it to values close to e and thereby 
satisfying the necessary conditions for BEC. 

One can give examples of phenomena from everyday life in which the BEC of quasipar- 
ticles in essence takes place. For example, from the standpoint of the quantum theory of 
solids, such processes as the ringing of various vessels, the sounding of tuning forks, etc., 
that are just the BEC of phonons "pumped" by the striking impact or, in other words, 
the preparation of a phonon Bose condensate, which, in turn, generates coherent sound 
emission. Remarkably, all this occurs at normal temperatures and so does not require 
any special experimental contrivances. At the same time, serious efforts continue to be 
made in the investigation of the BEC of such quasiparticles as excitons and biexcitons in 
semiconductors [9] and magnons in certain classes of quantum magnets [10]. Perhaps the 
most impressive are the recent results of Demokritov and co-authors [11] with microfilms of 
yttrium iron garnet. Their Mandelstam-Brillouin scattering experiments cleared revealed 
a resonance peak in the spectral density of the distribution of magnons over states as a 
function of frequency in the vicinity of the minimum energy of the corresponding energy 
spectrum upon an increase of pumping below a certain threshold value. These experiments 
might be described as having created an analog of the "magnon tuning fork", and at high 
(room) temperatures. 

Although it can be said that our present understanding of the various properties of 
Bose condensates of particles and quasiparticles is adequate, a number of unanswered 
questions remain, in particular, in regard to fluctuations of some observable quantities. 
This pertains primarily to systems of finite size or with a finite number of condensing 
particles, when the very concept of thermodynamic limit becomes problematical. Here it 
should be mentioned that the samples used in experiments on BEC are not only finite but 
usually small in size. If the spatial dimensions of the system are finite, then the statistical 
mechanics of the Bose gas is substantially complicated. In particular, the concept of 
equivalence of canonical ensembles loses meaning. We note here that in the BEC regime, 
inequivalence of ensembles is manifested even in the thermodynamic limit. Some finite-size 
effects for a Bose gas in a box were discussed in [12]. The "poor" behavior of fluctuations 
of the condensed particles described in the framework of the grand canonical ensemble 
(GCE) prompted the authors of [12] to draw the radical conclusion that this ensemble 
is unsuitable for describing any real physical system undergoing BEC. We add that the 
entropy also behaves "poorly" in the GCE - as the temperature goes to zero it does 
not go to zero in accordance with the Nernst theorem but, on the contrary, diverges 
(logarithmically) when the number of particles N — > oo . However, the GCE, because of 
its simplicity, is extremely convenient for concrete calculations; moreover, the behavior of 
such quantities as, e.g., the number of particles in the ground state, N , or the specific 



heat does not differ in calculations using different ensembles. 

The experimental study of BEC in traps has revived the theoretical research on 
different aspects of the statistical mechanics of Bose systems [13]- [23]. This is primarily 
because of the fact that until then attention had mainly been devoted to a spatially 
homogeneous gas found in a certain volume. In the traps that were actually used the gas 
is inhomogeneous, and therefore the results obtained previously for a homogeneous Bose 
gas in a box can be reproduced for a trap by considering it to be a potential well with a 
harmonic law of spatial confinement, i.e., basically reducing the problem to a system of 
N oscillators. 

Their partition function in the GCE (the grand partition function) is trivial to calculate. 
In the canonical (CE) and microcanonical (MCE) ensembles the corresponding partition 
functions are expressed in terms of contour integrals of the grand partition function. 
Therefore for not too large a number of particles and not very high energy, if one is 
talking about the MCE, the partition function of the Bose system can be calculated to 
any desired accuracy by the residue theorem or numerical integration. As to the analytical 
calculations of the partition function in the CE and MCE, one usually uses the saddle- 
point (steepest descent) method. However, as was shown in [12], for example, that method 
is inapplicable in the most interesting region - the neighborhood of the BEC, where the 
corrections to the main contribution do not fall off with increasing number of particles in 
the system. 

A number of calculations have been done to investigate the thermodynamic properties 
of finite Bose systems in different statistical ensembles. In particular, in [13] the difference 
between the behavior of the number of condensate particles N Q and also their fluctuations, 
calculated for a harmonic trap with the use of the GCE or CE when the total number of 
particles varied in the range 10 2 < N < 10 6 . In [14], which was devoted to a comparison 
of the results of an exact calculation in the MCE write approximate results obtained by 
the saddle-point method, it was noted that these results differ substantially precisely in 
the neighborhood of the BEC point. 

Nevertheless, to this day there is no convenient analytical representation for the 
partition function in the CE and MCE which would correspond well enough to the BEC 
regime. The goal of the present study is to remedy this. In Sec. I we introduce the necessary 
notation and definitions and also obtain an analytical expression for the average number of 
particles No on the ground level under the condition N ^> 1 . In the Sec. II we calculate 
the partition function in the CE by the saddle-point method with the first correction 
taken into account. In Sec. Ill we analyze why the domain of applicability of the saddle- 
point method is limited to temperatures T > Tbec '■> here by isolating the singularity 
corresponding to the ground level, we derive an expression valid for T < Tbec as well- In 
Sec. IV we propose a method of asymptotic expansion of the partition function in the CE 
in inverse powers of the number of particles, which works both above and below the BEC 
temperature. On the basis of the representation obtained, we calculate the fluctuations of 
the Bose condensate and demonstrate the quantitative agreement with the exact result 
down to very small values of the number of particles in the system. In the Conclusion we 
discuss the BEC temperature for a harmonic trap and a box and also the difference of 
the mathematical mechanism of formation of the critical point in the GCE and CE. 



1 Grand Canonical Ensemble 



A stationary quantum system consisting of N noninteracting particles is known to be 
completely characterized by the configuration [n] = {n , ni, n 2 , . . .} , where n^ = 0, 1, 2, . . . 
is the number of particles in the kth quantum states ((k = 0,1,2,...) According to 
the precepts of statistical mechanics, the (time) average of an observable quantity A 
in a nonstationary system coincides with to the average over an ensemble of stationary 
systems. An ensemble is determined by the distribution function p[n] . Then 

A = Z- 1 ^p[n]A[n], (1.1) 

[n] 

where the normalizing coefficient Z (the partition function) has the form 

Z = 5>H- (1.2) 

w 

In the description of an ideal Bose gas one generally uses the GCE, with a distribution 
function 

£n fc ( M - £fc )/T 

p[n] = e* , (1.3) 

where T is the temperature, fi is the chemical potential, and Sk is the single-particle 
energy of the fcth state. Here and below we have set Boltzmann's constant ks = 1 . Since 
p[n] (11.31) is factorized with respect to the dependence on the occupation numbers n k , 
the summation over configurations [n] is trivial to do. In particular, the average number 
of particles in the k th state 

oo 

ne n(/x-£fc)/T 

n fc = £ p[nW = ^ = e{£ J )/T _ 1 • (1.4) 



n=0 

The function on the right-hand side of ( 11.41) specifies the average occupation number and, 
hence, is a constituent element of the expressions for the majority of thermodynamic 
quantities (and not only in the GCE), and it is deviate to use a special notation for it: 

1 _ 1 

nk ~ e (e*-/*)/T _ l ~ e ( efc - eo )/T( 1 + _ l ■ tl-5) 

With the use of (5) the partition function (II. 2p is expressed as the product (II .5p 

oo 

Z=U(n k + l). (1.6) 

k=0 

The independent variables in the GCE are assumed to be T and \i . However, as follows 
from Eq. (11.51) . they could be considered to be the temperature and the number of particles 
n in the (ground) state with the lowest energy eo , which facilitates the analysis of the 
different regimes of the GCE. 

We write the average value of the total number of particles in the form 

oo 

N = J2^k = n + N ex {n T), (1.7) 

k=0 



where 



oo 



N ex (n ,T) = N cx = J2n k (1.8) 



k=l 



is the average number of particles in the excited states. If the value of N is fixed, then 
all the numbers n k except n fall with decreasing temperature (n^o = at T = 0), 
and no grows to tiq = N at T = . We denote by n k the maximum possible value of 
the average occupation number at a given temperature, i.e., the value of n k of Eq. (11.51) 
for n — > oo (or, equivalently, \i = Eq ): 

Then for no ^ 1 for the number n k specified by Eq. (11.51) one can limit consideration to 
the expansion 

n k (n k + 1) 

n fc ~ n k , (1.10) 

n 

which, in turn, reduces Eq. (I1.7P for n to a simple quadratic equation: 

_ ~ 5N 2 

N = n + N ex 2L, (1.11) 

n 

in which 

oo oo 

N ex = ^n" fe) 5iV e 2 x = ^n fc (n fe + 1) (1.12) 
fc=i fe=i 

is the maximum possible number of particles in excited states and its mean-square fluctu- 
ation. We note that the quantities marked with a tilde are functions of temperature only. 
The solution of equations (11.111) with respect to n is denoted as 



iVo(T) = I (^V - iV ex + ^(N-N ex y + ASNlj , (1.13) 

which we shall call the Bose condensate. This terminology is conditional in the sense that 
one is considering a problem outside the thermodynamic limit, with a finite total number 
of particles 77. It follows from the definitions ([TH]) and (fTT2l that N ex (T) and 5iV e 2 x (T) 
are monotonically increasing functions of temperature. We denote by T c the temperature 
at which A^ ex is equal to iV , which corresponds to 

N ex (T c ) = N. (1.14) 

If 5iV 2 x < (N - N cx ) 2 , and this holds for T ^ T c and N > 1 , then the behavior 
of solution (11.131) in the limit iV — >• oo acquires a stepped character. For different 
temperature regions, both below T c , where N > N ex , and above T c , where iV < N ex , 
the asymptotic behavior of N (T) at large but finite N has the simple form 

T<T C , 

N (T) = { 6N ex , _ T = T C , . (1.15) 

-iV), T>T C 




The value of T c at which the change of regime (the crossover) in the behavior of N (T) 
occurs can be regarded as a generalization of the temperature Tbec to the case of a finite 
number of particles in the system. We recall that Eq. (11.111) is approximate, in accordance 
with the condition no 1 . In the opposite case, when n < 1 (the Boltzmann limit), 
Eq. (|1.7I) gives the simple depends 

No{T) ~ 1 + Q(T) ' 3( T ) = E e " (£fc " £ ° )/T ' (1-16) 

which attest to the classical behavior of the Bose systems under consideration. 

Independently of the number of particles n in the condensate the factorized character 
of the distribution function ( jl.3f> in the GCE is conditional upon the absence of any 
correlations between particles of the Bose gas in different quantum states. This has the 
consequence 

n^n7 = nVn;. (1.17) 

The average of the square (and higher powers) of the number of particles in the k th state 
is calculated in analogy with Eq. dl .41) : 

rf k = 2n 2 k + n k , (1.18) 

from which the mean-square deviation (or, in other words, the mean-square fluctuation) 
is easily calculated and has the form 

5n 2 k = nf - n k 2 = n k {n k + 1). (1.19) 

Taking Eq. ( 11.171) into account, we write the square of the fluctuation of the total number 
of particles as 

oo 

6N 2 = Sn 2 + 6N 2 X , 6n 2 = n (n + 1), 6N 2 X = ^ n k (n k + 1). (1.20) 

k=l 

For T < T c the value n ~ iV , and the square of the fluctuation of the number of 
condensed particles is 

5n 2 = n (n + l) ~iV 2 . (1.21) 

Thus the relative fluctuation (Suq/ ' N) 1 / 2 grows with increasing number of particles in the 
system, and this is the basis for the widespread assertion that the fluctuations diverge 
below the BEC point (see, e.g., [2]). 

We note in this regard that the description of the BEC in the framework of the GCE 
cannot be considered quite correct, if for no other reason that it explicitly violates the 
Nernst theorem. Indeed, the entropy of the GCE is expressed in terms of the average 
occupation number as 

S = J2^ Hk + X ) ln ( Uk + 1 ">~ n k^n k }. (1.22) 

k 

In the region T<T C , where n ^> 1 and n k ^ <C 1 , it becomes equal to the entropy of 
the Bose condensate: 

S* ~ (n + 1) ln(% + 1) — no hin ~ ln(n + 1) + 1. 

When T — > , the entropy 5* ~ In iV , i.e., not only does it not go to zero but it diverges 
with increasing number iV . As will be seen below, in the canonical ensemble there is no 
problem with a divergence of the fluctuations nor with the entropy. 



2 Canonical Ensemble 



The main difference between the CE and GCE is that the total number of particles in the 

oo 

CE is rigidly fixed: N = ^ n fc = n + N ex . From this it follows directly that 

k=0 

n = N-N cx , (2.1) 
(n -n ) 2 = (iV cx -Ar ex ) 2 , (2.2) 

i.e., the fluctuation of the Bose condensate does not differ from the fluctuation of the total 
number of particles in excited states, or 

5n = 5N ex . (2.3) 

It is essential here that in the CE the average number of particles in the k th state is not 
equal to the average occupation number, determined in Eq. (11.51 ). and the noninteracting 
particles in different quantum states (in contrast to the GCE) are correlated with each 
other. 

The fluctuation of the Bose condensate in the CE can be estimated starting from the 
following qualitative arguments. For T > T c the number n <iV. Therefore, considering 
the Bose condensate as a small subsystem, one can suppose that a description of it in 
the framework of the GCE is valid. Consequently, for T > T c the following relation also 
holds in the CE [see Eq. (fTTTOl) ]: 

cm 2 , ~ n {n + 1). (2.4) 

When T < T c , however, the small subsystem becomes the particles above the condensate, 
and now the mean-square fluctuation of their number in the CE can be described by the 
expression 

oo 

5iV e 2 x ^n fc (n fc + l). (2.5) 

k=l 

Doing a simple interpolation of expressions (|2.4K and (|2.5p with (|2.3p taken into account, 
we find that 

oo 

n (n + 1) n k(n k + 1) 

S< = = . (2.6) 

£ n k{n k + 1) 

fc=0 

Below we compare this phenomenological expression with the expression calculated directly 
in the CE. 

For this we define the distribution function in the CE [cf. Eq. (11 .3ft 

p[n]=e * 6(N-J2*k). (2.7) 

k 

The property of factorization of p[n] is lost because of the presence of the 6 -function on 
the right-hand side of Eq. (12.71) . so that in comparison with the GCE the calculation of the 
partition function and the averages of the observables is complicated. Factorizability can 
be easily restored, however, and the summation over configurations [n] can be reduced to 



a summation over independent n^ if one uses the integral representation of the Kronecker 
5 -function: 

7T 

5(m) = — I dxe ixm . (2.8) 
2tt J 

—TV 

Then, substituting (12. 8ft into (12.71) and changing the sequence of summation and integration, 
we arrive at the following representation for the partition function: 



l - Idxe^Te-^^ = ±- hxe^ w ^ , (2.9) 
In J z — ' 271 J 



2tt 
where 

W(-ix) ln (! - e~ £k/T ~ ix ), 

k 

or, denoting x as iv , 

X 

W(u) = J2Hn k + l), n k = _ v _ l . (2.10) 

A comparison of Eqs. f )2. lOH and p. 61) shows that the function W^v) is the logarithm of 
the partition function in the GCE if one sets v = fi/T . Thus in the CE the averages of 
physical quantities are just the ratios of the corresponding integrals of the same averages 
found in the GCE. This, on the one hand, establishes a definite relation between the 
two statistical ensembles and, on the other, allows one to use the initial [see Eq. ( 11.51 )] 
definition for the number of Bose particles in each of the states. However, in them the 
average in the CE can also be expressed in terms of derivatives of the partition function. 
For example, in the simplest case of a nondegenerate spectrum e k we have 

n k = 5>[n]n fc = -Z~*T± $>M = -T^M , (2.11) 



de k 



Sn 2 k = -T^ . (2.12) 



At large the integral on the right-hand side of Eq. (12.91) can be evaluated by the 
saddle-point method, which leads to the following asymptotic expansion: 



e W(i>)-vN 



v) 



'l + Zl + z 2 + •••)■ (2-13) 



Here v denotes the saddle point nearest to the origin of coordinates in the complex 
x -plane (x s = iv) , the equation for which has the form 

oo 

W'{u)=^2n k = N. (2.14) 

fc=0 

In the leading asymptotic approximation the logarithm of the partition function in the 
CE has the following simple form: 



\nZ = W{v) -vN-\ \n[2itW"(u)\. 



(2.15) 



The contributions Zj in Eq. fj2 . 13f ) are expressed in terms of ratios of the derivatives of 
the function W{y) of the type 

[W®{y)r[WW{y)]» (2 lfi x 

[W"(v)] m+n+ i ' ' ' ' 

If the function W(u) and its derivatives are large, W®(v) ~ iV , and this is the 

least in the region T > T c and N 3> 1 , then the ratios ( 12.161) and the contributions Zj 

have order of smallness 0(N~i) . Thus for the first correction we find: 

_ WW(V) 5[W"\u)r _ 

We take the derivative of the occupation number f12.10f) with respect to Si : 

T^ = n k (n k + l)(r^-5 k ^. (2.18) 

The derivative of the saddle point v with respect to ei is evaluated by differentiating 
Eq. ((214)): 

m dv nAni + 1) . 

Now with the aid of Eqs. (j27TTI> . (127T5D . (I2T8D . and (12191) we find the average value of 
the number of particles in the k th state, 

n k (n k + l) ( W"(u)\ 

E ' = ^-^W( 2 "' +1 -Fiy)- (2 ' 20) 

In particular, for the average number of particles in the ground state (12.111) we obtain 

E ° = n ° ~ n fw»t? [{2no + 1)V>) " V "' {v)] ' (2 - 21) 

where V{v) denotes the sum over only the excited states, 

oo 

V{u) = W{v) - ln(n + 1) = ^ ln(n fe + 1). (2.22) 

k=l 

finally, differentiating n with respect to e [see Eq. (|2.12l) ]. we arrive at an expression 
for the mean-square fluctuation of the Bose condensate: 

6n 2 = 6 1 + 6 2 + 6 3 , (2.23) 

with the leading contributions 

52 = 2W\vY { {2n ° + 1)V "' {V) ~ (6n ° + 6n ° + 1)V " {V) + (2 - 25) 
+ 2 ^^ +1) [(2no + l)V>)-^)]}, 



2V"'(u] 
W"(v) 



[(2n + i)y»-y'"(v)]}. 



In the temperature region T > T c the condensate is dilute, n <C N , and, accordingly, 



Si ~ n (n + 1), 5 2 




(2.27) 



We note that the phenomenological formula (12.61) for the fluctuation of the Bose condensate 
coincides with the leading asymptotic contribution 5i (12.241) . 

An exact expression for the partition function is given by the single integral (12.91) . 
which for not too large N is easily found numerically. It is interesting here to compare 
the exact expression for the fluctuations with its asymptotic behavior given by formulas 
(I2.23I) - (I2.26I) . Such a comparison, however, is impossible to do in general form, since the 
quantitative calculations require specifying the explicit form of the function W(u) , which, 
in turn, depends on the concrete form of the energy spectrum Sk ■ Let us find it for the 
case of alkali metal atoms in magnetic traps. 

Experiments on cooling of a large number of alkali metal atoms (N ~ 10 3 . . . 10 4 ) are 
interpreted as the experimental realization of BEC. The particles are confined in the traps 
by a potential v(r) , the exact dependence of which on the distance r is, gives speaking, 
unknown, but for theoretical analysis usually a quadratic (harmonic) approximation is 
used. As a result, the problem of BEC reduces, as we have said, to a calculation of the 
partition function of a system of linear oscillators. The spectrum E\ and the spectral 
density g\ (coefficient of degeneracy) of the three-dimensional isotropic oscillator has the 
simple form 



We point out that in Eq. (12 .281) the index I enumerates the energy levels and not quantum 
states, which are enumerated by the index k introduced previously. At high temperatures 
(T e ) the series expressions for the function W{v) and its derivatives converge slowly. 
It is shown in the Appendix how to improve their convergence and to obtain expressions 
convenient for numerical calculations. 

The equation for the crossover point in the CE differs from the equation (11.141) in the 
GCE by only the replacement of A by A , i.e., 



It follows from definitions (TTTT21 and (l27[4l) that A ex = W'(f) - n . Then, using for 
iy'(z/) the asymptotic expansion [Eq. ( 1A.21I) in the Appendix] for high temperatures 
r = T jhuj ^> 1 , we write the following expansion for A ex : 



where £(j) is the Riemann ( -function. From it we can find the solution of equation (12 . 29f) 
that determines the crossover point r c in the form of an expansion in inverse powers of 
N . In the leading approximation we denote this solution as 




(2.28) 



A ex (T c ) = A. 



(2.29) 



A ex = r 3 C(3) + \r 2 C,(2) + r lnr + O(r), 



(2.30) 




N 



T 1/3 



TBEC = 



C(3) 



(2.31) 



Now, knowing equation ( 12.311) . we find for r c from (12.291) and (12. 30ft 



r BEC 



BEC 



C(3)r BE c 



C(2) Iiitbec 
2 3r B Ec 



0(t, 



BEc)' 



(2.32) 



It is seen from expressions (12.311 ) and ( 12.321) that the crossover temperature T c is 
below the condensation temperature Tbec for a Bose gas in a trap. We note that for a 
Bose gas in a box the situation is the opposite, T c > Tbec • 

Expressions (12.311) and (I2.32H with the known numerical values of the ( -function in 
them easily convince one that even for a number of particles of the order of 10 3 the 
condensation temperature T c is only 6 times greater than the ground-state energy e . 
This is a direct indication that the condensation phenomenon observed in the experiments 
mentioned is of a microscopic (or, in any case, mesoscopic) rather than macroscopic 
character. This casts doubt on whether the condensation of several thousand particles 
can be regarded unambiguously as BEC, the main feature of which, strictly speaking, is 
the appearance and manifestation of quantum properties in macroscopic phenomena or 
objects. 

Without denying, of course, the presence of the phenomenon of BEC itself in magnetic 
traps, we would nevertheless like to say that, in our view, the results set forth in this 
Section are evidence that the condensation of alkali metal atoms observed in the experiments 
is more of a nanophysical character. 

In Fig. 1 we show graphs of the fluctuations of the Bose condensate in the case of their 
exact calculation, 

^ d 2 InZ 



5nl 



del 



and their approximate calculation ( 12 . 231) — (12 . 261) . It is seen that the asymptotic expressions 
represented by Eq. ( 12.231) and the corrections to it do not adequately reproduce the curve 
of the numerical calculation. 



3 Bose-Einstein Condensation Region T < T c 

Let us consider in more detail the low-temperature region, where one can more or less 
definitely talk about the presence of a Bose condensate. As we have said (see Fig. 1), 
here the discrepancy between the exact and asymptotic values of the fluctuations are 
significant. The reason is not hard to understand. The fact is that in the region T < T c the 
evaluation of integral ( 12.91) by the straightforward saddle-point method does not actually 
give the asymptotic expansion in inverse powers of N: the terms Z\ , Z2, ■ ■ ■ in Eq. ( 12.131) 
do not fall with increasing N . The contributions of the ground term Wo to the sum 

oo 

W(v) = Wk (12.10P and its derivatives with respect to v are: 

w = ln(n + l), w' = n , w'q = n (n + 1), w'q = n (n + l)(2n + 1). 

Since for T < T c the pth order derivative ~ ~ N p , for the first correction, 
e.g., we have 

_ l w (rio + l)(6rig + 6wo + l) _ 5 ng(n + l) 2 (2n + l) 2 _ _J_ 
Zl "8 n 2 (n + iy 24 n 3 (n + l) 3 ~ 12' 




Phc. 1: The dependence on temperature t = T/Tbec of the relative fluctuation of the Bose 
condensate Siiq/N . The circlets denote the results of a numerical calculation of the integral in 
Eq. (|2.9p . Curve 1 is the leading contribution to the asymptotic expansion (I2,23p . curve 2 is with 
the next correction to Eq. (12.231) taken into account, and curve 3 is the asymptotic expansion 
HTTP . (I4T41) . 



It is not hard to show that the other contributions Zj in (12.131) not only do not fall with 
increasing N but even grow with increasing j . Nevertheless, the problem of singular 
behavior of the contributions corresponding to the ground state can be solved as follows. 
We denote by U(—ix) the exponent of the integrand in Eq. (12.91) : 

U(-ix) = ixN + W(-ix). (3.1) 

It follows from the definition of W{y) [see Eq. (12.101 )] that the function U(—ix) is singular 
at the points xi — i£i/T , and its real part goes to infinity at these points. Of course, the 
function U(—ix) reaches its minimum value at points iv\ lying along the imaginary axis 
between each pair of singular points: 

u'(ui) = o, e ^-<vx< e ^. 

The behavior of the function U{y) in the vicinity of the first singularity xo is shown 
schematically in Fig. 2. 
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Phc. 2: Behavior of the function U(u) in the vicinity of v = Eq/T . 



The depth of the minimum of the function U{y) at the saddle points vq and V\ 
depends on temperature: U{v Q ) < U(ui) for T > T c and, oppositely, U(uq) > U(i>\) for 
T < T c . Therefore for an optimal estimate of the integral <\2M for T < T c the integration 
contour must be deformed so that it passes through the saddle point x s = iv\ , as shown 
in Fig. 3. At the point x the function exp[U(—ix)] has a simple pole. The contribution 
from this pole to the integral (I2.9p . which we denote Z , is equal to the residue of the 
integrand there: 



where the divergence V(u) is defined in Eq. (12.221) . The contribution of the integral along 
the contour C has the form 



the contributions from the parts of the integration contour along the ray [— 7r, — n + zoo) 
and [it, n + zoo) , being equal in magnitude and opposite in sign, cancel each other out. 
Furthermore, it is easy to see from expressions (13.21) and (|3.3I) that for T <T C and N ^> 1 
the ratio of the integrals Z\jZ§ is exponentially small, and therefore the partition function 
in the BEC regime is given by the exceedingly simple expression (13.21) . It follows from 
that expression, in particular, that the entropy 
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We note that, as a consequence of the periodicity of the function W{y) , i.e., 

W(v + 2m) = W{v), 
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Phc. 3: Integration contour for the integral in Eq. (|2.9p for T <T C . 



goes to zero at T — > , as it should. 



4 Modified Asymptotic Expansion 
of the Partition Function 

The most interesting region, but the hardest for calculations, is the critical neighborhood 
of T c . Here the contributions of Z and Z\ are of the same order of magnitude, and 
the fluctuations are maximal. To obtain the correct asymptotic expansion of the integral 
(12.91) in inverse powers of the number of particles N we propose the following approach, 
consisting of several steps. 

i) In the first step we separate out explicitly the first singular term in the integrand 

7T 

(4.1) 
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where, in contrast to (13.11) . the function U(-ix) here has a different form: 

U(-ix) = | + a(iV+i)+y(-a). (4.2) 

The saddle point x s = iu for the function (14.21) satisfies the equation [cf. Eq. (12.141) 

V'{y) = N + \. (4.3) 
ii) In the second step we make the change of integration variable x — > u : 

u 2 = u(u) - U(-ix) = V(u) - V(-ix) - (u + ix)(N + i), (4.4) 



ludu = iU'(-ix) = i[V'(-ix) - N - \\. 



(4.5) 



Hi) Finally, we deform the integration contour in the x -plane so that it passes through 
the saddle point along the line of steepest descent, which is determined by the equation 

lm[U(-ix)} = 0. 

As a result of these steps we can transform the integral in (14. ip to a form in which 

Z = Z + Zt, (4.6) 

oo 

1 (4.7) 
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Zt = B — J du e- u 'f(u), (4.8) 

— oo 

and denote by iv the value of the variable u , defined in Eq. (14.41) . corresponding to the 
point x = ieo/T : 

v 2 = U(f) - U{v) = V(f) - V{v) -(f-v)(N+ \). (4.9) 

The function f(u) in the integrand of (14.81) is analytic in the neighborhood of the 
point u = iv: 

/(") = -rm ",..m m 1 Z^ITaT ~ 77T~-77 ' (4-10) 
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where the function x(u) in fj4. lOj) is determined by Eq. (14.41) . The integral in (14.71) can 
be expressed in terms of the error function: 

Z = \ e^+^Erfc (v) = \ e v{£o/T) - N£o/T Eik (v), (4.11) 

which satisfies the well-known equations [24] 

oo 

Erfc (v) = I dx e~ x \ Erfc (0) = 1, Erfc (-oo) = 2, Erfc (-v) = 2- Erfc (v), 
V 71 " J 

(4.12) 

and, for v 3> 1 , it has the asymptotic expansion 

Erfc ( v ) = ±^ (l - J- + JL + 0(t;- 6 )) . (4.13) 
Vv/vrV 2f 2 4f 4 7 



We find an approximate value of the integral ( 14.81) by expanding the function f ]4. 1 OH 
in a Taylor series at the point u = : 

^ = ^[/(0) + i/ ff (0) + •••]■ (4.14) 

It follows from Eq. (14.51) that dx/du = —2iu/U'(ix) . Ultimately we obtain for the function 
f{u) and its second derivative at the point u = 

1 1 
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In the limit T — > T c we have v — > Eq/T , v — > 0. Then, resolving the uncertainty 
(oo — oo) at the point v — , we obtain from (14.151) for the crossover region 
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The contributions /(0) , f"(0) , and their ratio /"(0)//(0) all fall off with increasing N . 
The power of the decrease of the ratio /"(0)//(0) depends on the asymptotic behavior 
of the derivatives V^(u) , which, in turn, is determined by the concrete form of the 
spectrum E\ or the spectral density gi as functions of I . For the quadratic confining 
potential considered in the Appendix, the derivatives ~ N for p < 3 , ~ A^ln N 
and ~ A^ P//3 for p > 3 . It follows from this that representations ( 14.111 ) and f!4. 14H for 
the contributions to the partition function are indeed asymptotic expansions in inverse 
powers of N for all temperatures, including the critical temperature region T ~ T c . 

As we see from Figs. 2 and 4, even at relatively small numbers of particles (as low as 
N ~ 100 ) the discrepancy between the exact expressions for the fluctuations of the Bose 
condensate and the approximate expressions corresponding to the representation of the 
partition function in Eqs. (14.61) , ( 14. lip . ( 14 .141) , and (14.151) is indiscernable on the graphs. 
It can be shown that at a temperature below T c the given representation coincides with 
the asymptotic expression (13.21) . and for T > T c it goes over to fj2 . 13H . 

The critical region is determined by the condition 

|Erfc(u) -1| ^e, e<l. (4.17) 

For example, setting e = 10 _1 , we find from ( 14.171) the boundary value v e = 1.16 and, 
using the asymptotic approximation for V(v) [see Eq. (1A.211) in the Appendix], we obtain 
for the critical temperature region \T — T c \ ^ AT 
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(4.18) 



It is seen from Eq. (14. 1 8j) that this region narrows extremely slowly with increasing particle 
number. In other words, the thermodynamic limit is reached at very large particle numbers 
iV > 10 6 . Figure 5 illustrates the evolution of the relative fluctuation of the density of 
the Bose condensate as the total number of particles in the system is varied. The curves 
shown convincingly demonstrate that even for a system with iV w 10 4 the temperature 
behavior of these fluctuations nevertheless differs markedly from the limit iV — > oo , in 
which 

fag V"(e /T) n _t(2)( T \ 3 

5 Conclusion 



As we have said, in the standard theoretical treatment of the BEC problem the trap is a 
three-dimensional box with, importantly, a fixed volume. In this case the Bose condensate 




Phc. 4: The difference between the approximate expressions for the fluctuations and the exact 
values (<5ng exact — <5nQ as )/iV. Curve 1 corresponds [see Eq. (|4.15l) ] to the leading asymptotic 
contribution /(0) , curve 2 includes the asymptotic correction /"(0) . 



formation temperature (on the quantum scale r = T / ftuj^^) has the form 

r ^-(JL.Y' 3 , 51 ) 

while in the problem considered above, for a trap in which the confining potential is 
quadratic and the volume is not fixed, 

1/3 



« = ( 7^ ) • (5-2) 



N 

As is seen from Eqs. (15. II) and (15 .2fi . the dependence of the BEC temperature on the 
number of particles in the system is significantly different in these two cases. For the box 
the characteristic energy is expressed in terms of the volume V of the system and the 
mass M of the particle in the following way: 

n 2 h 2 , s 

^ box = 2^7^ • (5 " 3) 

In the magnetic trap the volume is not fixed and, moreover, it changes with changing 
temperature. We define the effective size of the system in a trap with a quadratic confining 
potential as the amplitude of the oscillations of an oscillator with energy equal to the 
temperature T . Then 

V^-(— Y 7 *. (5 4) 
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Phc. 5: Dependence of the relative fluctuation 5uq/N on temperature t = T/Tbec for different 
values of the number of particles in the system. 



Now, after expressing the number of particles in terms of the particle density p (N = p V) 
and substituting Eq. (15.31) into (15.ip and (15.41) into (15.2ft . we find for T BEC in the two 
cases 

2*2 / \ 2/3 
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(5.5) 
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It follows from (15.51) that in the language of particle number density the the temperatures 
of BEC in the box and trap are close not only qualitatively but also quantitatively: 



T B b EC _ 1 p^C(3) 
T^l 2 [4C(3/2) 
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-1.13. (5.6) 



On the whole, it should be emphasized one again that BEC is rightfully considered to 
be one of the fundamental discoveries of theoretical physics. Its clearest trait is not simply 
the accumulation of a macroscopic number of particles in the ground state upon cooling 
of an ideal Bose gas but the fact that this process is a phase transformation in a system 
of mutually noninteracting particles. Most likely the term "Bose condensation" came into 
use because of the analogy (which, strictly speaking, is not entirely correct) with the 
condensation of a vapor to a liquid, which was discussed by Einstein in his pioneering 
paper [3]. 



Although phase transformations in nature are extremely diverse, at the same time 
they demonstrate surprising universality: the change of the thermodynamic properties of 
a system occurs abruptly when the temperature (or some other controllable parameter) 
crosses its critical value. From a formal theoretical physics point of view the main question 
is, how, in functions which are initially analytic in temperature, does the singularity arise 
at the critical point T c : 

f(T,N) = h{T)6{T c -T) + f 2 {T)e{T c -T). (5.7) 

N— ioo 

Bose-Einstein condensation, as an exactly solvable model, gives a simple answer to this 
question. For example, in the GCE the role of the "smeared" 6 -function, according to Eq. 
(I1.13p . is played by the quantity 

Mi) ^( vx- + L/ ')' —-S- (5 - 8) 

and in the CE, as follows from Eq. (I4.1ip . by 



9 CK {x) = ^Erfc (-a CA VNx), a CA = J (5.9) 

For both representations ( 15.81) and ( 15.91 ) at iV — > oo the limit is the ordinary 6 -function, 
but at a finite value of N they behave differently. It is seen in Fig. 6 that ^ca is closer 
to a step than 9 GC a ■ Thus one can say that the thermodynamic limit sets in somewhat 
faster in the CE than in the GCE. 

If as the function f(T, N) in (15.7P we take the specific free energy 

f(T,N) = ~hiZ, 

then both fi(T) and f2{T) are identical in the two ensembles, thus implying the equivalence 
of the ensembles: 

/i(T)= "c(3)^' /2<T)= "^(3r^- (5 - 10) 

The chemical potential /i in Eq. (15.1 OH as a function of temperature and particle number 
density is the solution of the equation 

T 3 Li 3 (e^)=C(3)T| EC , 

where Li k (z) is the polylogarithm [see Appendix, Eq. (1A.12|) ]. However, the specific fluctu- 
ations of the Bose condensate is different in the different ensembles: in the GCE 

l im ^1 = / 00 nP H T < T c 

n^oo N \ npn T > T c , 

and in the CE 
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Phc. 6: Behavior of the smoothed # -functions for N = 10 2 . Curve 1 — @ca{x) , curve 2 — 
0gca{x) . 



In and of itself the accumulation of particles at the lowest energy level is the direct 
and rather trivial consequence of Bose statistics, evident from the form of the formulas for 
the average occupation numbers, Eq. (11.51) . When one is talking about a phase transition, 
however, the question of thermodynamic limit must also be addressed: for example, is a 
number of particles iV = 10 3 sufficient to approach it? To speak more precisely, the term 
thermodynamic limit is commonly understood to mean letting the volume of the system 
go to infinity at fixed temperature with the various densities (e.g., the particle number 
density p) held constant. However, the volume, density, and temperature are dimensional 
quantities: 1 meter is almost "infinite"if one is measuring in angstroms. >From the point of 
view of theoretical physics, any limit should be formulated in the language of dimensionless 
quantities, and in the problem of BEC of an ideal gas there are only two such quantities: 
the total number of particles in the system, N , and the temperature on the quantum 
scale, r = T/ftxj. Therefore the thermodynamic limit is, first and foremost, N — > oo , 
and the densities should be taken as the ratios of the corresponding extensive quantities 
to the total number of particles. The asymptotic corrections to the limit N — > oo in 
the our problem are of the order of iV -1 / 3 , so that, it would seem, the thermodynamic 
limit is reached with 10% accuracy if N — 10 3 . On the other hand, however, the BEC 
temperature for this number of particles is comparable to the ground state energy e (we 
recall that in the given problem T c ~ 6eo at N = 10 3 ), and so one cannot speak of a 
macroscopic scale of the physical quantities. 

In this regard we note that when atoms of alkali metals are held in magnetic traps 
the procedure of preparing a coherent state of N particles is said to involve "cooling" 
apparently to reflect its thermodynamic nature. Meanwhile, neither the volume nor the 
temperature nor, moreover, the spectral density of the particle number cannot be controlled 



to the required precision because of technical shortcomings and for fundamental reasons: 
what is the temperature on quantum scales T ~ £ ? Therefore, in light of the results 
presented above, there is ample justification for concluding that the assertion that true 
BEC has been observed in these undeniably outstanding experiments should be taken 
with a degree of caution. 
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Appendix 

The function W(v) (12.101) and its derivatives W^'iy) is determined by series of the form 

oo 

^2 h{l) . Separating m the first m terms, we denote the remainder of the series as 



1=0 



l=m+l 



and take an Abel-Plana transform of it. Then the series in Eq. (I A. II) is transformed to 
the sum of two integrals, 

K = I + J, (A.2) 



in which 



dlh{l), (A.3) 
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J = "2 L k-^ -■ (A.4) 







In order for series ( 1A.1I) to converge, the function h(l) must fall off with increasing / 
not slower than / _1 . For a power-law function h(l) ~ l~ a its derivatives hS p \l) behaves 
at large / as l~ a ~P . In this case the integral J can be evaluated with the aid of an 
asymptotic expansion in inverse powers of m . We expand the function h(m + | — ix) 
in the integrand of ( 1A.4I ) in a Taylor series at the point x — . For the sum of any 
finite number n of terms in this expansion one can switch order of the summation and 
integration and get 



n 

J~Y(-l) l c l h { - 2l + 1 \m+ A), (A.5) 
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where 



1 f dxx 2l+1 (l-2- 2 ^ 1 )C[2/ + 2] 
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We note that the optimal number of terms in the asymptotic expansion f lA.51) is determined 
by the form of the coefficients c\ . Therefore evaluation of the integral J can be done with 
the aid of (1A.5I ) to any desired accuracy by increasing m . 

We now use this technique to evaluate the logarithm of the partition function (12.21) . 
Here the function h{l) has the form 

h(l)=g(l)Hni + l), (A.7) 

ff(/) = ^(l + !)(* + 2), m = e(l _J /T _ 1 , r=£, o = vr-\. (A.8) 

Then 

W» = ^g(l)ln{n l + l) + I + J. (A.9) 

The integral I in fl A.9[ ) is expressed in terms of the polylogarithms (Lerch functions) 

/ = r 3 Li 4 (^) + ar 2 U 3 (z) + brU 2 (z), (A.10) 

where 

a = g'(m + \) = m + 2, b = g(m + \) = \(m + |)(m + §), 

z = e -(-+|--)A ; (A.H) 
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The contribution J in (1A.9I) has the form 

n 

J=J2(-l) l c l f l , (A.13) 
«=o 

where 

/o = aln(r + l)-6r -1 r, (A.14) 

for I > 

/, = -Z(2Z + l)r- M+1 r«_i + a(2Z + l)r~ 2 V 2i - br^W, (A.15) 

1 t d l r ^ dr i-i 

r = r = n m+ i = — — - L — , n = t — = r(r + 1)— — . 

For calculating the asymptotic expansion of W{u) at large r it is sufficient to keep 
only one term in the sum on the right-hand side of (IA.9I) : (m = in (1A.1I) ). Then, using 
the well-known asymptotic expansion of the polylogarithm for x — > 

Ii4(e-) = C(4) - C(3)x + C(2)^ + Anx - ^ ^ + 0(x 4 ), (A.16) 
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where is the Riemann (^-function 



C(4) = ^, C(3)^ 1.202, C(2) = y 



with accuracy up to terms that do not fall off with increasing r we find that 

I = ^ + g'(a)asy + g(a)^T+ 



90 

1/2-Q- 



{ [g(a + 1) + §] ln[i^] - § [ll(<r + 2) 2 - a - § ] } + 0(r^] 
Here the contribution (1A.4I ) takes the form 
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/dx i 
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(A.17) 
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The asymptotic expansion for the derivatives W^\v) we obtain by differentiating expres- 
sions (1A.10I) and (IA.18I) with respect a : 

, n T P f) P 

W^(v) = —[Hn + l) + I + j}. 

We note that, starting with the third derivative, expression (IA.10I) for / can be written 
in terms of elementary functions, since 

Lii(e- x ) = -ln(l-e- x ). (A.19) 

The derivative dJ/da is expressed in terms of the special function ip(z) = T'(z)/T{z) : 

1 



00 r- / 1 

dJ „ f dv r ; ' 



-2 / — Im 

da J e 2wx + 1 
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I — a — ix 



|s 0(a) [V(l - <x) - ln(| - a)] . (A.20) 



In particular, for the derivative W'{y) = 9(l) n i with Eq. (|A.20l) taken into account, 

1=0 

we arrive at the expression 



Ti 2 



W\v) = n + C(3y + g , (a)-r 2 + Tg(a)[\nr-ij(l-a)] + 

+l(3a 2 + 5a - f ) + 1(1 - a) [g(a + 1) + f] - 1 + O^ 1 ). (A.21) 
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